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Quantum information processing using quasiclassical electromagnetic interactions 

between qubits and electrical resonators 
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Electrical resonators are widely used in quantum information processing with any qubits that are 
manipulated via electromagnetic interactions. In nearly all examples to date they are engineered to 
interact with qubits via real or virtual exchange of (typically microwave) photons, and the resonator 
must therefore have both a high quality factor and strong quantum fluctuations, corresponding to 
the strong-coupling limit of cavity QED. Although great strides in the control of quantum infor- 
mation have been made using this so-called "circuit QED" architecture, it also comes with some 
characteristic limitations. In this paper, we discuss a new paradigm for coupling qubits electromag- 
netically via resonators, in which the qubits do not exchange photons with the resonator, but instead 
where the qubits exert quasi-classical, effective "forces" on it. We show how this type of interaction 
is similar to that induced between the internal state of a trapped atomic ion and its center-of-mass 
motion by the photon recoil momentum, and that the resulting multiqubit entangling operations are 
insensitive both to the state of the resonator and to its quality factor. The method we describe is po- 
tentially applicable to a variety of qubit modalities, including superconducting and semiconducting 
solid-state qubits, trapped molecular ions, and possibly even electron spins in solids. 

PACS numbers: 



Microwave electrical resonators are already an impor- 
tant tool for manipulating quantum information using 
electromagnetic interactions [l[ . As a means of quantum 
information storage or communication, they are used or 
proposed in a variety of schemes involving, for exam- 
ple, trapped molecular ions neutral polar molecules 
[3j], Rydberg atoms 0, Q, superconducting Josephson- 
junction qubits 0, 0] , and electron spins in solids [1, Q . 
In all of these cases, the resonator is used in a way famil- 
iar from optical cavity quantum electrodynamics [l(| , 
in which the qubits exchange real or virtual photons with 
the cavity, and where the figure of merit is the speed with 
which this exchange occurs (or equivalently, the strength 
of the coherent qubit-cavity coupling). In fact, many of 
the seminal cavity QED results have been replicated in 
a circuit environment, a field now known as circuit QED 
0,13 (cQED). 

Although this paradigm has been extremely fruitful for 
quantum information processing, it also presents some 
characteristic limitations in this context. For example, 
photon exchange can in most cases be controlled dynam- 
ically only by detuning the qubit from resonance with 
the cavity [11], which does not in general allow a very 
strong suppression when the coupling is intended to be 
off. In addition, whenever exchange of photons is possible 
between the qubit and other circuit elements, including 
but not limited to resonators, additional pathways open 
for the qubit 's excited state to decay, by emitting pho- 
tons into environmental degrees of freedom associated 
with those elements [l2[. Finally, since cQED exploits 
a strong, manifestly quantum interaction between qubit 
and resonator, and effectively uses the resonator to store 
or transport quantum information, its protocols are sen- 
sitive both to spurious population of excited resonator 
quantum levels, and to spontaneous decay of those lev- 



els. High fidelity in cQED quantum operations therefore 
requires both a correspondingly high cavity quality factor 
Q, as well as strong suppression of the spurious popula- 
tion of excited levels which is typically encountered in 
experiments [l3j . 

In this paper, we describe a different paradigm for 
using electrical resonators to couple (possibly distant) 
qubits via electromagnetic interactions, which does not 
involve any photon exchange between resonator and 
qubits, and therefore largely avoids these limitations. 
Our proposal is built on an analogy with demonstrated 
methods 0, SI developed by M0lmer and S0rensen [l6[ 
and Milburn [l7| for entangling trapped atomic ions via 
their collective center-of-mass vibrational modes [l8|, [lj| . 
In contrast to cQED as described above, here we manip- 
ulate only quasi-classical resonator states whose vacuum 
fluctuations can be small, using purely classical, qubit- 
state-dependent forces [20j. As a result, the resonator 
can be treated classically, its initial state is almost ir- 
relevant to the protocol, and high quality factors are no 
longer required. 

In section [I] below, we describe the general qubit- 
resonator system under consideration. Section [IT] con- 
tains a detailed description of how entangling operations 
can be achieved between qubits coupled to a common 
resonator, without any exchange of photons as in cQED. 
In section IIIII we consider in detail the leading sources 
of error in two-qubit entangling operations, and evalu- 
ate these errors for a variety of different physical qubit 
modalities, which are tabulated in table |U We conclude 
in section IIVI with a comparison to the methods used 
with trapped ions, as well as cQED. In the appendix we 
describe the assumptions and parameter values used to 
evaluate the gate errors shown in the table. 




FIG. 1: Electromagnetic qubit / resonator coupling without exchange of photons, (a) shows the electric circuit analyzed in detail; 
(b) is its exact dual [22I |23| . governed by equations of identical form, via the transformation: Qo$,C«i,VttI,ytt2 
[24| . in which loop charges [red arrows in (a)] become node fluxes [red dots in (b)] [2^|. The qubit states |e) and \g) have 
stationary energies E q (p) and E e (p), classical functions of p £ Qq,& q , and are described in (a) by a parallel state-dependent 
voltage source V q and capacitance C q (defined in the text), and in (b) by a series state-dependent current source I q and 
inductance L q . Resonator damping, and the associated fluctuations, are modeled in (a) via R r and the Langevin noise source 
6Vr, and in (b) by G r and SI r . (c) and (d) show classical systems analogous to the circuits of (a) and (b), in which an electric 
or magnetic dipole is placed inside the resonator, producing an orientation-dependent displacement, (e)-(j) show examples 
of physical qubits coupled to resonators in the manner of our proposal (the qubit bias is not shown) as well as the energies 
E g (p) and E e (p) schematically [see the appendix for details], (e) tunable transmon 0,0 (f) tunable flux qubit [26j, H3] (g) 
singlet-triplet double quantum dot [2^ - 4301 ] : (h) polar molecular ion where the third curve (shown in magenta) corresponds 
to the two degenerate J = 1, rtij = ±1 levels; (i) 31 P coupled electron and nuclear spins in 28 Si [3lfl. where we show only one 
nuclear spin orientation. The splitting at zero field is due to hyperfine coupling; (j) Nitrogen- vacancy center in diamond [321 ] . 
with a fixed transverse magnetic field, as a function of an additional longitudinal field. The third (magenta) curve shows the 
two degenerate levels which become ms — in the limit of large longitudinal field. 



I. GENERAL QUBIT- RESONATOR SYSTEM 

We first consider a microwave resonator coupled to a 
single qubit, as illustrated in figs.[TJa) and (b) for capac- 
itive and inductive qubit/resonator couplin g. T hese two 
circuits are chosen to be exactly dual [2% |23|. so that 
the solution for one case can be mapped directly to the 
other using the transformation: Q -H> <&, C f-> L, V O 
J,yHZ. The qubit-resonator interaction of interest to 
us is of a quite general form, in which the qubit's station- 
ary energy levels \g) and |e) have energies E g {p), E e (p) 
with a specified, classical dependence on the parameter 
p, which is typically the flux through a qubit loop or 
the charge on a gate electrode, induced by the resonator. 
We are interested here in cases where any interaction by 
which the qubit can exchange photons with the resonator 
can be neglected, several examples of which are shown in 
figs. [He)-(j). We consider the capacitively-coupled cir- 
cuit of fig. [Ha), and write the classical potential energy 
in terms of the loop charges Qb,Qc,Qr [25| (for the mo- 
ment taking R r , SV r — > 0): 



U(Q r ,Qc,Qb) =E q (0) 
V q (Q' q )dQ' q 
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V b (Q' b )dQ' b 



(1) 



where E q (Q q ) is the qubit energy, Q q = Qb + Qc, and 
V q (Q q ) = dE q /dQ q is the voltage across it; Vb(Qb) is the 
bias voltage which corresponds to Qb- The terms in the 
second line are the potential energies of the qubit and 
the source, respectively, and can be viewed as work done 
by (or on) the source as Vb is turned up from zero. 

We now seek to write eq. Q] in terms of deviations of 
the charges Q n (n € b, c, r) from the minimum energy 
solution Q n0 , given by: Q b a = C b {V b - V q ), Q c q = 
—C c V q , Q r o — Qcq (which must be obtained self- 
consistently since V q depends implicitly on Q q ). Note 
that U(Q r o,Qco,Qbo) is a constant, independent of V b , 
as required by total energy conservation. Writing: Q n = 
Qno + &Qn, we then expand the qubit energy: 



E q {Q qQ + SQ q ) ^ E q {Q 



qO) 



(2) 



where C q = 



E q (Q q )-E q (0) 



(d 2 E q /dQ q )^ 1 is the effective dynamic ca- 
pacitance of the qubit at Q q o- We emphasize that V q 
and C q depend on the internal state of the qubit, in 
addition to its bias. Of the three SQ n , only 5Q r is an 
independent dynamical variable, since it couples to an 
inductance (the effective "mass" for a fictitious parti- 
cle whose "position" is 5Q r [H, [25T ? ]), while SQb and 
6Q C are deterministically related to 5Q r . To determine 
these relations, we hold Vb and 5Q r fixed and minimize 
U with respect to SQ C and SQb, to obtain: 5Qb = and 
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SQ r /C r = SQdC- 1 + C~ x + C- 1 ). Combining this with 
eqs. [T][5J we hnd: 



= (Or + C c V q ) 2 



with the qubit-state-dependent quantities (via C q [33j): 



k', which imparts a recoil momentum fi(kx — k') (radia- 
tion pressure) and can be used to implement dissipativc 
laser cooling of the atomic motion [34j; or (ii) via co- 
herent (stimulated) scattering far from resonance (the 
so-called "dipole force" [35|), which can be used to real- 
ize essentially non-dissipative forces with internal-state- 
dependence. Here we focus on the latter case |36[. 



C' = C r 



C c Cq 

C C + Cq 



\J L r C' r 




T~ = ^T (4) 



Thus, the qubit effectively exerts a state-dependent 
"force" on the oscillator, displacing its equilibrium "po- 
sition" by an amount C c V q , without changing the total 
energy at that position. It may at first seem strange 
that a quasistatic displacement of Q r is possible at all, 
since L r would appear to act as a short across C r at low 
frequencies (for R r ,SV r = 0). However, Q r (which we 
choose to be the canonical position variable of the res- 
onator), is not simply the charge on the capacitor C r . 
Rather, Q r is more precisely a quasicharge [23J, defined 
in this case by: dQ r /dt = I r + C r dV r /dt, whose conju- 
gate momentum is the inductor fluxoid $ r defined here 
by: d& r /dt = V r — L r dl r /dt. A nonzero quasistatic dis- 
placement of Q r is therefore possible due to a nonzero 
time-integral of the displacement current C r dV r /dt. 

To see intuitively how this arises, we can map our prob- 
lem onto a more familiar electrostatic system, as shown in 
fig-[H c ) : that of an electric displacement C c V q contained 
inside the capacitor C r . We can imagine this as an elec- 
tric dipole which was moved into the capacitor (with its 
orientation held fixed), such that the resulting transient 
voltage appearing across C r gives J C r {dV r /dt)dt = C c V q 
(analogous arguments apply to the magnetic case in fig. 
QJb) for a qubit with dE/d& = I q ^ 0, which acts like 
a magnetization inside the resonator inductor L r , as in 
fig. Did)). This mechanical analogy casts the work done 
by the source (the second line of eq. [IJ into the form of 
a potential energy associated with "assembling" the sys- 
tem, as is often encountered in electro- or magnetostatic 
problems: as a permanent electric dipole is moved into 
the resonator capacitor, there will be a mechanical force 
on that dipole due to the gradient of potential energy in 
real space, such that work must be done while moving 
the dipole into position. In our system, this "motion" 
consists of turning up the bias voltage Vb- 

The form of eq. [3] can now be used as the basis for 
an analogy between our coupled qubit / resonator system 
and a harmonically trapped atom whose internal states 
are coupled to its center-of-mass motion by the photon 
recoil momentum. In the latter case, this coupling occurs 
typically in one of two ways: (i) via Rayleigh scatter- 
ing of laser photons tuned near resonance with an elec- 
tronic transition from an incident laser wavevector kf, to 



II. CONTROLLED-PHASE GATE WITH 
QUASICLASSICAL FORCES 

Following refs. we describe the oscillator in 

terms of the dimensionless classical field amplitude: a = 
(it + iv), where: 



Q r a + fit 




(5) 



and $ r is the resonator flux (associated with L r ), the 
canonical momentum of the oscillator such that [u, v] — i. 
The quantity a is a c-number for coherent states, while in 
the presence of statistical (e.g. thermal) fluctuations of 
the field it can be written as a diagonal density matrix (in 
the Glauber-Sudarshan P-representation [37[). When we 
couple N qubits to our resonator, in general all resonator 
quantities become operators in the 2 Ar -dimensional qubit 
Hilbert space: a — > a, u> r — > t2> r , Z r — > Z r , Y r — > Y r . In 
the cases of interest to us here, these operators are all 
diagonal, and commute with the individual qubit Hamil- 
tonians [381 ] . 

We wish to focus on the qubit-state-dependent dynam- 
ics of the operator a, while the state dependence of the 
other resonator quantities are higher-order effects which 
we will consider as perturbations only, in section Hill To 
that end, we define: 



LU r = Q r (l + Su> r ) 



(6) 



where ui r = Tr[p m ui r ], and p m = I 2 n /2 is the 
completely mixed state of the N qubits (Id is the d- 
dimensional identity matrix). This definition separates 
the average effect on the resonator due to coupling to the 
qubits (which is simply a renormalization of its frequency 
and impedance) from qubit-state-dependent effects. Wc 
can now write the equation of motion for the oscillator 
with quality factor Q as [391 ]: 



dot ct — ck* 

_ = + (7) 



where r = uj r t, we have made the following replacements 
in eqs.[5]above: Z r = 1/Y r — > Z r = l/Y r = Lo r L r , and we 



can write the dimensionless, 
quency shift and force as: 



qubit-state-dependent fre- 



N 



6w r = ^ 5uj~&* 

i=l 

AT 

fj = [(T + &rr)6i + (V + + Sr, + %] (8) 



where i € 1...N indexes the qubits, of and Ii arc the 
Pauli-z and identity operator for qubit i, and: 



C' 



C C + Cn 



C C V:: 



(9) 



using the notation: X ± = [(e|X|e) ± (g\X\g)]/2 with 
X a single-qubit operator. For the sake of clarity we 
have assumed all N qubits and coupling elements are 
identical (though the method we propose does not require 
this), and we have retained only the leading order term in 
Sw~ <C 1 (we will see in table U below that this is a very 
good approximation) . Note that there is no nonzero <5w+ 
because of our definition of ui r [c.f., eq.[S]. The quantities 
Si) are Langevin noise terms, with Sr)~ due to qubit bias 
noise, and Srj + associated with the finite resonator Q. In 
our model [fig. QJa)] the latter comes from the Johnson- 
Nyquist noise SV r of the resistance R r (Q = Z r /R r ). 
The resulting dimensionless noise power spectral density 
of the fluctuating force 5i) + can be written: 



20, 



'-(SV r (t)SV r (t')) u 
Q,T r . 



Q 



coth 



(10) 



where the brackets denote an environment average, the 
subscript u) indicates a Fourier transform, O = io/uj r is 
dimensionless frequency, and r c = hiu r /kBT r with T r the 
effective resonator mode temperature. 

To implement a multiqubit gate, following the trapped 
ion case [T3 . [l6| , we use time- and qubit-state-dependent 
forces on the oscillator, which result in qubit-state- 
dependent paths of the oscillator in u, v phase space. 
Each amplitude in the A^-qubit wavefunction then ac- 
quires a geometric phase associated with the correspond- 
ing oscillator path, given by the phase space area it en- 
closes: <j) g = lm[§ a* da]. To produce the forces, we 
drive all of the qubits simultaneously with an oscillatory 
bias [13]: Qbo — Qdc + SQAC s ' m {^dt), such that the 
corresponding state-dependent force on the oscillator is: 




FIG. 2: Oscillator phase space trajectories due to classical 
"forces" exerted by two qubits. Trajectories are shown in a 
frame rotating at ui r , and the area enclosed within each path, 
indicated by solid lines (\gg) - blue, \eg),\ge) - black, |ee) - 
red) is the geometric phase acquired by that amplitude of 
the qubits' internal-state wavefunction. Dashed arrows show 
the displacements of \ee) and \gg) relative to \ge), \eg), to 
make a connection between the general case considered here 
where (eg\f)\eg), (ge\fj\ge) ^ 0, and the usual situation in the 
trapped-ion case {eg\fj\eg), (ge\fj\ge) = [T3 - [l^ ]. (a) In the 
ideal case (Q — > oo, 8u) r = 0), the paths are closed. The in- 
set shows the high-frequency component of the response due 
to the counter-rotating term in eg. 1111 (b) for finite Q,8ui r , 
the trajectories are no longer closed, such that at the end of 
the gate a nonzero total oscillator displacement 5ctq is en- 
tangled with the qubits' internal state. The inset shows the 
components 5u^f , Sv~ of this residual displacement (5v£ = 
because 8ui^ = 0); finite decay results in a shrinking radius 
of the path with time, producing finite 8u~ at the end of the 
gate, and finite 8u) r makes the drive detuning 8d weakly qubit- 
state-dependent, producing a nonzero Svq . (c) modified gate 
sequence where either the force is reversed, or the qubits are 
inverted, in the middle of the gate, such that at the end of the 
gate the spurious entanglement between qubits and oscillator 
is removed. 



fj = fjo sm(u)dt) (for a drive centered on a bias point where 
V~ = 0, a so-called "degeneracy point" [41]). Neglect- 
ing for the moment the finite Q, fluctuations Sr) , and 
frequency shift Sui~ , and taking the drive to be turned 
on at r = with the oscillator in the state cto and not 
entangled with the qubit, we find (in a frame rotating at 



a = iNUoe lSdT 
f)o 



iS d T 



s d 



5 d (2 



(11) 



where the detuning 6# of the drive is defined according 
to: aid = £v(l + Sd). In eq. II H the first term is free evo- 
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lution, and the second and third terms can be viewed as 
the co- and counter-rotating components of the oscillator 
response, respectively, due to the effective force f] ac - 

The resulting oscillator dynamics in phase space, for 
a two-qubit system, are shown in fig. Ufa). The entan- 
glement of the oscillator with the qubits vanishes when 
the trajectories return to their initial point, at times r 
satisfying the conditions: 



2tt 



T 

2^ 



(12) 



where m, k are integers {k > 2m). During this time the 
system has accrued a state-dependent geometric phase, 
given by (neglecting overall phases): 



2nm 



s d 



A controlled-7r gate, sufficient in combination with single- 
qubit gates for universal quantum operations, is imple- 
mented (in addition to the single-qubit phase shifts given 
by the second term in eq. [T5]) if we choose: 



2 + 8 d 
4m 



(14) 



where the total gate time is, to leading order in S d : 



2tt 



m 1 
2 Vo 



(15) 



Obviously, choosing m = 1 gives the fastest possible gate 
for a given 77^; however, as we will see, this choice in 
most cases does not give the lowest gate error. 

Before discussing realistic errors in these operations, 
we highlight some distinctions between the type of cou- 
pling we are discussing, and that used in cQED @, Q- 
First, the coupling strength [c.f., eq. [5] here depends 
only on the resonator admittance, and not on its reso- 
nant frequency, as compared to cQED where the detun- 
ing of the qubit from the resonator controls the strength 
of the coupling. This difference is a direct consequence 
of the fact that in our scheme no photons are exchanged 
between resonator and qubit. In addition, here the cou- 
pling strength is inversely proportional to the vacuum 
fluctuation amplitude of the oscillator y/HY r , such that 
the larger the classical displacement of the oscillator com- 
pared to its quantum fluctuations, the stronger the cou- 
pling; By contrast, in circuit QED, the opposite is true: 
for example, in the case of a transmon qubit coupled ca- 
pacitively to a resonator, the transverse coupling strength 
can be written: g cx (3 r \/hY r - strong interaction corre- 
sponds to the case where quantum fluctuations dominate 
and classical displacements are absent (this is precisely 
the strong coupling limit in which circuit QED becomes 



possible) . An important consequence of this is that in cir- 
cuit QED spurious population of resonator excited states 
(equivalently, nonzero resonator temperature) linearly af- 
fects gate fidelity, whereas in our scheme it has only a 
higher-order effect for statistical fluctuations, and no ef- 
fect for coherent resonator excitation. This difference 
may be important since not only is it difficult to com- 
pletely remove spurious excited-state populations [l3j . 
but it is also much harder to achieve very high resonator 
Q at the single-photon level than it is when the resonator 
is strongly driven [42| . 



III. TWO-QUBIT GATE FIDELITY 

We now seek to evaluate the constraints on the two- 
qubit gate fidelity, and examples of our results are tab- 
ulated in table Q]for several qubit modalities. We begin 
with the effect illustrated in fig. [U[b), spurious entan- 
glement between qubits and resonator at the end of the 
gate in the form of a residual qubit-state-dependent os- 
cillator displacement Sa~ . The two mechanisms which 
produce this type of entanglement are finite resonator Q, 
and the state-dependent frequency shift Su> r . As shown in 
fig-HUb), resonator decay causes the radius of the phase- 
space trajectory to shrink with time, while 8ui r effectively 
produces a detuning Sd (and corresponding evolution rate 
around the paths in phase space) which is qubit-state- 
dependent. We can bound the resulting error eg a using 
the fidelity [43| between the desired qubit density matrix 
after the gate pf and the actual total density matrix 
after the gate pj?*, traced over resonator states jiij ]: 



£<5q 



tot _ n t 
Pf ~ Q 



Do = D 



Tr 



yJipfyWrlpTKP?) 1 ' 2 



pf ® D\a a + 5a+)p r th D(a + 6a+) 



(16) 
Dq 



\gg)(gg\+D 



5a„ 



|ee)(ee| 



where p r th is a thermal resonator state with tempera- 
ture T r , D(Sa) is the displacement operator in u, v phase 
space for (complex) displacement 8a, and Dq performs 
the state dependent displacement shown in fig. [2](b). 
Note that in the presence of damping, «o is nonzero only 
if the resonator is driven. As an approximate worst case 
estimate, we take: pf — (\gg) + \ee))((gg\ + (ee|)/2, and 
find (using the results of ref. \4m : 



e<5< 



\5a q 



-(l + 2n) 



(17) 



to leading order in 5a~ and n ~ exp (— hu> r /kBT r ), the 
mean thermal photon number in the resonator, and in- 
dependent of a + 8o>+ . 

Now, for the gate shown in fig. [U(b), the finite Q 
and 5Cj t result in displacements: Su~ ~ 0(l/Q) and 
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5v~ ~ 0(5uj~), respectively. However, since both of 
these displacements result from classical and determin- 
istic dynamics, we can strongly suppress them using the 
scheme shown in fig. Oc). We divide the gate excitation 
into two periods of equal duration (each of which satis- 
fies eqs. [T2"l) . and we switch the sign of the effective force 
rj ac between these two periods (alternatively, one could 
insert 7r-pulses to invert the qubits instead). This sim- 
ple procedure works like a spin-echo, in the sense that 
the classical modifications to the trajectories cancel out 
at the end of the operation, such that the leading terms 
~ 0(1/Q), 0(5u>~) cancel. To evaluate the result quan- 
titatively, we first renormalize the oscillator resonance 
frequency to account for the usual Q- induced shift: 



1 

4Q2 



(18) 



and correspondingly renormalize the dimensionless time 
r [eq. [7], drive frequency uJd, detuning Sd, and the con- 
ditions for m and k [eq. [12]. The final result, to leading 
order in the small quantities 5d, T)q , Su>~ , is: 



\5a 



-\ 2 



1 + 2m 



(19) 



where x = ^/{QfJo) and rn now includes both halves of 
the gate. As we show in table HI this echo procedure re- 
duces es a to very low levels for all of the cases considered. 

In addition to the displacement Sv q , the state- 
dependent frequency shift 8u~ also produces direct de- 
coherence due to thermal photon-number fluctuations. 
Since this shift can be viewed as an n-dependent qubit 
frequency: uj q b + 2nSuj~ , each qubit experiences dephas- 
ing at the rate @: 



sy 16ui r nQ(Su> r Y 



(20) 



for a thermal photon number distribution at temperature 
T r with mean photon number n [46| . Notice that the 
dephasing rate increases with increasing resonator Q, as 
the discrete frequencies associated with different photon 
number states become more resolved @, [l3| . 

We now consider the classical fluctuations Srf 1 of the 
force during the gate, which arise both from the damping 
and possibly from qubit bias noise. For small fluctuations 
about the desired values, the system will accrue an addi- 
tional geometric phase: 

6<f> g « 2 J h V - [ V + (ai + a z 2 ) + »T*"*a] 

+ Sr) + [r)-(af + a z 2 )]\dT (21) 

where the three terms are: single- and two-qubit phase 
errors due to fluctuations of the state-dependent force 



(qubit bias noise), and single-qubit errors due to oscilla- 
tor fluctuations (as well as qubit bias noise if V q + ^ 0). 
Since rf 11 oscillate at aid, and only noise which occurs at 
or near this high frequency or its harmonics will pro- 
duce errors, in nearly all cases low-frequency qubit bias 
noise (e.g. 1/f charge or flux noise) can be ignored in 
this context. The oscillator noise, however, can be quite 
important as the Q is reduced, producing an average er- 
ror (per qubit): 1 - F « (6<%) = / S Sr ,+ {Q)\Ti-(Q)\ 2 dn 
where r/~(£l) is the fourier transform of tj~(t), a peaked 
function centered on fid, of width ~ (2t/) and ampli- 
tude ~ r lo T f- We obtain: 



(22) 



These errors will tend to restrict how small a Q can be 
used for these gates [47| . 



For qubits without a degeneracy point, such as singlet- 
triplet quantum dots |28l - t3Qj | , there will be an additional 
error contribution from the Q or $ sensitivity associated 
with a nonzero rj, in the presence of noise in these param- 
eters. Of course, this type of error will always be present 
for such qubits independently of our scheme, so that de- 
coupling schemes will be required (48j in any case. Here, 
the simplest such scheme would be to replace the reversal 
of r\ in the middle of the controlled- 7r gate with a 7r-pulse, 
exchanging the roles of \g) and |e) for each qubit (sim- 
ilar to ref. 30). In the limit where the 7r-pulse is much 
shorter than the gate time, this will produce the same 
effect as shown in fig. [U[c), as well as the usual spin-echo 
to remove the effect of quasistatic noise [13, EH • 

Looking back at eqs. nTll2"2l we conclude that: (i) larger 
forces 77^" are better; and (ii) higher-frequency resonators 
are better, for faster gates, lower n, and so that hui r 
ksT (r c 3> 1). Some specific examples are shown in 
table |H where we attempt to illustrate the utility of 
our scheme over a range of qubit types, resonator Q's, 
resonance frequencies, and photon populations n. Note 
that in some cases two-qubit error rates as low as 10~ 3 
are still achievable with Q as low as 25,000 or n as high as 
0.4, illustrating the robustness of our technique against 
resonator decay and fluctuations. With state-of-the-art 
resonator Q's of 10 6 and low n, even lower error rates ~ 
10~ 4 become accessible (provided that single-qubit errors 
unrelated to the operations considered here do not limit 
the gate). 

In this context it is important to note that while very 
high resonator Qs of 10 6 or greater are at present difficult 
to achieve in planar geometries at the single-photon level 
required for cQED implementations, this difficulty is sub- 
stantially reduced at higher resonator powers (~ 10 3 — 10 4 
photons) where parasitic, lossy defects become saturated 
|42| . Given that our proposed method is not affected 
by a coherent initial resonator state, we can consider the 
prospect of intentionally driving the resonator to increase 
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qubit 


Q 


uj r /2n 




Cr 


V 


8 d 


m 


ri 1 

<? 


^77 






£2qb 


Aim 


iV 7 


type 




[GHz] 


r tti 


[P F ] 


[x 10 3 ] 


[xlO 3 ] 




r l 

M 


r 1 

[ns] 


[xlO 3 ] 


[x 10 3 ] 


[x 10 3 ] 


[GHz] 




Flux 


25000 


10 


250 


0.84 


4.5 


64 


100 


8.4 


160 


0.015 


2.0 


2.0 


2.3 


2600 


[26] 


10 6 


10 


500 


0.46 


3.4 


12 


6 


0.76 


51 


10" 9 


0.27 


0.34 


3.5 


5200 




10 6 


1 


30nH 


0.84 


1.4 


2.9 


2 


0.23 


700 


2 x 10" 8 


2.0 


5.0 


0.03 


3.1 x 10 6 


Transmon 


50000 


10 


250 


0.84 


1.4 


14 


50 


126 


360 


0.053 


4.5 


4.6 


0.4 


10 4 


[12] 


10 6 


10 


250 


0.84 


1.4 


6.2 


10 


6.3 


160 


7 x 10" 8 


0.51 


0.53 


0.9 


10 4 




10 6 


3 


1500 


1.8 


1.1 


2.2 


2 


0.10 


300 


3 x 10~ 8 


1.5 


4.5 


0.07 


2.0 x 10 5 


S-T QDs 


25000 


15 


lOOnH" 


l.lfF 


3.0 


33 


60 


3.0 


120 


0.20 


3.8 


4.0 


1.4 


2.5 


[28-301 


10 6 


15 


30nH a 


3.7fF 


1.7 


6.6 


8 


0.83 


80 


5 x 10~ 7 


0.47 


0.57 


1.4 


8.5 


CaCl+[2j 


10 6 


1 


1000nH a 25fF 


0.14 


1.1 


30 


25s 


27^s 


0.0032 


5.1 


5.1 


0.01 


2.4 x 10 9 


31 P- 28 Si [31] 


10 7 


1 


100 b 


250 


0.014 


0.14 


50 


2 x 10 7 s 


350/is 


0.005 


4.0 


4.0 


0.5 MHz 


2.4 x 10 9 


NV [32] 


10 7 


1 


100 b 


250 


0.028 


0.13 


10 


2s 


78^ts 


6 x 10" 5 


4.4 


4.5 


0.8 MHz 2.4 x 10 9 



a Impedances this high require the use of high-kinetic-inductance materials [541I55H . Limits on achievable resonator Q and impedance 

for these materials are as yet unknown. 
b Impedances this low may not be achievable in a transmission-line resonator. 



TABLE I: Selected examples of two-qubit controlled- ir gate parameters (Details for each system are contained in the 
appendix). For transmission-line resonators, L r and CV are effective values for a given longitudinal mode. The results for 
gate errors shown are the thermal photon- number dephasing rate F^ [eq.[20] (which is present due to the coupling even when 
no two-qubit gates are being driven), the gate time t^ [eq. [15], the state-dependent displacement error es a [eqs. 117119] and 
the resonator fluctuation error e Sv + [eq. [22]. The total two-qubit error is defined as: t2 q b = 2r t f,t^+es a +^s v + ■ This does not 
include single-qubit errors unrelated to the coupling (for example, Ti relaxation for the superconducting qubits or charge- 
noise dephasing for the quantum dots). The quantity Ahm is the minimum detuning to the next higher oscillator mode, if 
we assume it has the same Lamb-Dicke parameter as the fundamental, such that we the associated error ehm < 0.1 x e2 q b 
[eq. 123] . In all cases we take a resonator temperature of T r — 40 mK, corresponding to n — 0.4,0.03,4 x 10 _6 ,10 -8 
for uj r /2n = 1, 3, 10, 15 GHz. Note that both donor spins in Si and NV centers in diamond do not themselves require 
low temperatures, but here they are required to suppress errors due to classical, thermal resonator fluctutations [c.f., 
eqs.ilEU]). 

the effective Q for gate operations (note that the ther- 
mal fluctuation errors discussed above associated with T r 
would also apply to the effective temperature of an addi- 
tional drive field). The simplest way to accomplish this 
would be to drive the resonator instead of the qubits, 
and use the qubit /resonator coupling to accomplish the 
modulation of f) for the gate. In table Q] we list the quan- 
tity 7V 7 , the resonator photon number (at oj^) required to 
produce the full qubit bias swing necessary for each set 
of gate parameters. In nearly all of the cases considered, 
these values are comparable to or larger than than the 
~ 10 3 — 10 4 photons typically necessary to saturate the 
resonator loss in current experiments. 

One disadvantage of this scheme is that by driving the 
resonator, all qubits connected to it would necessarily 
become entangled. In architectures where this would be 
a problem, an alternative method would be to drive the 
resonator with a separate, resonant drive (i.e., at £o r ) to 
saturate its loss, while driving only the desired qubits 
at ujd', to avoid spurious entanglement from the resonant 
drive in conjunction with a nonzero 5uj~ , one would need 
to invert the sign of this resonant drive in the middle of 
the gate, resulting in a classical "echo" similar to the 
method discussed above and shown in fig. Etc). 

The final issue we need to address is the presence of 



higher modes in the resonator [12], which will also be 
displaced by interaction with the qubits, according to 
their effective capacitance (inductance for the magnetic 
case of fig-IUb)) an d geometrical coupling factors. The 
resulting (state-dependent) excursions do not in general 
decouple from the qubit at the same time as the funda- 
mental mode, leaving a spurious entanglement between 
the qubits and the higher modes. Working in our fa- 
vor, as in the atomic case is the fact that the 
excursion of higher modes is suppressed by detuning, by 
at least a factor 8d for transmission-line resonators with 
equidistant harmonics, and potentially significantly more 
for the lumped case. Using eq. [T7]for the error due to a 
state dependent excursion Sa k of mode k with frequency 
u! r< k = fk^r, and taking Sa k to be ~ n~ c k the effective 
Lamb-Dicke parameter for that mode [c.f., eq. ITT]: 

2 (l ™ fk ) 

Using this result, we list in table UJ the minimum sep- 
aration A^ m from the fundamental to the next higher 
mode, assuming that rj^ k = r/^ , which would give rise to 
an error ehm one tenth of the total error from all other 
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sources. In all cases A^ m < oj r , indicating these errors 
arc unlikely to be significant in practical cases. 



IV. CONCLUSION 

The method we have presented here has some favor- 
able advantages compared to the analogous techniques 
demonstrated with trapped atomic ions. For example, 
in our case, the resonator modes are nearly (though not 
completely) independent of the qubits, so that adding 
more qubits does not generate additional nearby col- 
lective modes which must be avoided as in the case 
of trapped ions. Also, in the ion case the equivalent 
controlled-phase interaction requires the ions to be in 
magnetic field sensitive states, which produces inevitable 
dephasing 14]; in our proposal, although the states used 
must also be field-sensitive, the sign of their sensitivity 
is oscillating at a high frequency so that to leading order 
dephasing is nearly absent (except in qubits with no de- 
generacy point). Although alternative methods for ions 
have been demonstrated that use states without field sen- 
sitivity [l5|, , they have additional complications which 
limit the gate speed, fidelity, and number of ions that 
can be entangled [1J|. Finally, although the schemes 
used for ions are nominally insensitive to the state of 
the resonator, they still require a small Lamb-Dicke pa- 
rameter, meaning that the ions must be localized to a 
much smaller region than the laser wavelength to avoid 
sampling a spatially dependent force [l6j]. This can be 
particularly challenging in the context of the ubiquitous 
heating observed in ion traps, which causes 77 to increase 
in time, in the absence of active laser cooling [U. In 
our case, although 77 is also a small parameter (in the 
expansion Eg^ e (p)) 1 this expansion does not break down 
as in the atomic case; in fact, most of the error sources 
we have discussed actually decrease with larger rj. In the 
expansion of eq. [5J we included terms up to second order 
in the resonator displacement 5Q ri and for the parame- 
ters in table Q]the errors associated with the second order 
term [c.f., eas. 119120] are already small enough for low er- 
ror rates. Under these conditions, third and higher-order 
terms have a negligible effect [50j . 

Our method also has some potentially useful ad- 
vantages over the current cQED paradigm for qubit- 
resonator coupling Since no photons are exchanged 

between qubits and resonator, there is no Purcell effect 
[l2[, and the qubits' excited-state lifetimes are decou- 
pled both from the qubit-resonator detuning, and from 
the resonator Q. It also avoids the relatively small on/off 
coupling ratios encountered in cQED architectures where 
the coup ling is controlled only by detuning qubits from 
cavities [ill l5lj . Furthermore, unlike in circuit QED 
where spurious population of excited photon states of 
the resonator directly produces gate errors, our scheme 
can tolerate even substantial thermal occupation of ex- 
cited resonator states (up to n = 0.4 was considered in 
table|T|, as well as coherent driving of the resonator, with 



low impact on gate fidelity in the cases we have consid- 
ered. This removes the need for an extremely low effec- 
tive resonator temperature and/or a high resonator fre- 
quency (i.e., T r fiWr/kB), and also in some cases the 
need to operate the resonator in the single-photon limit 
where high Q is more difficult to achieve. 

In closing, we note that the geometric phase of eq. [Rj] 
can be generalized to TV qubits [c.f., eq. 0, giving an ef- 
fective interaction cx (X^^f) 2 - This type of many-qubit 
interaction may be of interest, for example, in cluster- 
state generation [52| , or syndrome extraction in quantum 
error-correction schemes where multiqubit operators are 
required, such as surface codes 

Note added - An error-suppression method closely re- 
lated to that shown in fig. |2{c) was recently demonstrated 
for trapped atomic ions in [D. Hayes, S.M. Clark, S. Deb- 
nath, D. Hucul, I.V. Inlek, K.W. Lee, Q. Quraishi, and 
C. Monroe. Phys. Rev. Lett. 109, 020503 (2012)]. 

We acknowledge helpful discussions with and/or com- 
ments from Daniel Greenbaum, Arthur Kerman, and 
Adrian Lupascu. This work is sponsored by the Assis- 
tant Secretary of Defense for Research & Engineering 
under Air Force Contract #FA8721-05-C-0002. Opin- 
ions, interpretations, conclusions and recommendations 
are those of the author and are not necessarily endorsed 
by the United States Government. 



Appendix: Details on error estimates 

1. Superconducting qubits 

For the flux qubits [fig. QTf)], we take Ej/h = 200 
GHz, E c /h = 5.7 GHz, and a(<J> 60 ) = 0.74. This gives 
w gb ($ h0 )/27r = 4.73 GHz and L~ = 15.9 nH, L+ = [5^. 
We assume a drive amplitude of 8&ac = 0.1$o, corre- 
sponding to a(0.13> ) = 0.70 and cj g6 (0.1$ ) = 7.8 GHz. 
For the transmons [fig. QJe)], we take uj q b($bo = 0)/2ir = 
15 GHz, which gives L~ = L+/2 = -2($ /tt) 2 /hw p = 
—87.4 nH. Here we assume a drive excursion 5&ac = 
0.2$ , corresponding to w 9fc (0.2$ ) = 13.5 GHz. For 
both flux and transmon qubits we take L c — 25 pH [57] ]. 
Also for both of these cases we must consider the effect 
of junction asymmetry in the DC SQUID (defined by the 
area asymmetry parameter: A j = 2(Ai— A2)/(Ai+A2)), 
which produces a spurious coupling between external flux 
and the SQUID plasma mode. For the limit L c <C Lj 
satisfied here, we find the matrix element for the drive 
to couple to the first excited vibration al st ate of the 
plasma mode: Ma ~ (S^AC'/^o)Ajy / Tr/Ahajp. We 
then have a maximal probability in this excited state of 
~ (^/4)(<5$ AC /$ ) 2 ^(a; p /K - Lo d ))\ Thus, for Aj < 
0.05 and <5<J>ac/$o < 0.2, this error is at the 10~ 4 level 
or below. This calculation also gives a residual, direct 
Jaynes-Cummings type coupling between DC SQUID 
plasma mode and the resonator: g ss LOpfi r Aj ' Z r /8Rq, 
which can be neglected in the cases considered here. 
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FIG. 3: Orientation of vector qubits driven by rotating fields, 
(a) trapped molecular ion qubits, and (b) electron spin qubits 
driven by a rotating field Ed or Bd, respectively, whose rota- 
tion plane contains the resonator mode field direction; in both 
cases the dipole orientation follows the drive field. In (c), the 
NV center spin triplet is driven by a rotating magnetic field 
Bd whose rotation plane contains the resonator mode field 
(B r ) direction; the orientation (quantization axis) of the NV 
center is determined by the crystal, and is also contained in 
the drive field's rotation plane. 



2. Singlet-triplet coupled quantum dots 

Figure HJg): The singlet-triplet coupled quantum dots 
(and molecular ions below) are biased with a voltage, 
rather than a charge, so that: Cb,C c 3> C q (the op- 
posite limit from the superconducting qubit cases). In 
this limit, the voltages across Cb,C c can be neglected, 
and V q s» V&. The quasicharge is then: Q q ps dE~ /dV q 
evaluated at V/,. The displacement of the resonator in 
eq. |3] is then C q Vb and instead of eqs. [5] and [HI we 
have: % = {dE~ / 'dV q ) / \2^/W r ) , C~ = d 2 E-/dV q 2 , and 
5uj~ /bj r m —C~/(2C r ). For the quantum dots we use 
the parameters of ref. l28l with a tunneling amplitude of 
t c = 23/jcV. We take an exchange energy J for each qubit 
which oscillates from ~ 1/j.eV to ~ 2.5/_teV, corresponding 
to electrode voltages from -0.7mV to -0.4 mV, dE~ /dV e 
from 0.002e to O.Ole, and C~ = 2C+ -1 to 10 aF. 

3. Trapped molecular ions 

Figure dpi): The trapped polar molecular ions of ref. 
have a vector dipole moment, associated here primarily 
with the J = 1 excited molecular rotational manifold 
with sublevels m,j = — 1, 0, 1 which are degenerate at zero 
electric field. An electric field shifts the J = 1, mj = ±1 
levels down, so that a qubit can be realized with the J = 
0, mj — and J — 1, mj — states. We therefore cannot 
simply oscillate the electric field through zero, since the 
mj = state will undergo Majorana-like transitions to 
the mj — ±1 states near zero field (equivalently, the 
induced dipole's orientation will not "follow" the applied 
field) . Instead, we can use a rotating electric field Ed with 
angular frequency u>d, whose plane of rotation contains 
the resonator mode field axis as shown in fig-EJa). In this 
case, as long as the rotation is not too fast, the molecular 
dipole will follow it, resulting in an oscillating projection 
of the dipole along the resonator mode field. The effect of 
the rotation can be expressed via Larmor's theorem as an 
effective magnetic field along the rotation axis: B rot = 
hijJd/lJ=i where 7j = i is the gyromagnetic ratio of the 
,7=1 manifold. As long as the rotation is turned on and 



off slowly, and we restrict ourselves to ftwd ^ E m j= o — 
E mj= ±i, the states will transform adiabatically. 

With this in mind, we take Ed = 2.5 kV/cm, corre- 
sponding to a differential stark shift between the J = 
0, mj = and J = 1, mj — rotational states from their 
zero-field splitting of 9 GHz to ^13 GHz, and a splitting 
between J = 1, mj = and J = 1, mj — ±1 of ^3 GHz; 
we take a distance I — l/_tm between the qubits and drive 
electrodes, so that Vb ~ EJL. 

Finally, although the rotation at constant field magni- 
tude Ed does not couple the J = 0, mj = to the J = 1 
manifold, there is a residual direct Jaynes-Cummings 
type coupling between the resonator and the J = <-> 
J = 1 transition, given by: hg rj (doi/l)y/hY r /C r = 0.4 
MHz (dot ~ 2.5 x 10 _29 C-m is the transition dipole be- 
tween J = and J = 1); this produces a negligible effect 
at the large qubit-resonator detunings >10 GHz consid- 
ered here. 



4. Donor electron spins in Si 

Figure [Hi): Similar to the case of quantum dots and 
molecular ions which are biased with a voltage or elec- 
tric field rather than a charge, an electron spin is biased 
with a magnetic field rather than a flux. In this case, 
we have: % = (dE~ / dl q ) / \2^KZ~ r ) , L q = d 2 E-/dI 2 , 
and 5uj~ /uj r ~ —L~/(2L r ). We take a circular loop of 
diameter d — 50nm connected to the resonator, with the 
spin at its center, so that the field at the spin is given 
by: Bd — fiolq/d. As in the case of molecular ions, a 
rotating field (turned on and off slowly) must be used 
to avoid Majorana transitions between spin orientations 
[fig- [3jb)] , whose rotation frequency in this case should 
be less than the Larmor frequency = ^ e Bd/h. For 
our drive at 1 GHz, we then select a field amplitude of 
B d = 1000G, corresponding to u; l /2tt = 2.8 GHz. At 
this field, the hyperfine splitting of ~ 117 MHz [58| pro- 
duces a negligible L~ m 10~ 22 H. Our chosen parameters 
Q = 10 7 , d = 50 nm, and B d = 1000G, while arguably 
not completely implausible, are admittedly extreme. The 
very weak interaction with a single spin dictates that 
such parameters are required for favorable gate parame- 
ters. One way to relax these requirements to some extent 
would be to implement a multiturn coil (it would need 
to be at the ~ lOOnm scale) to increase the coupling, or 
to use an ensemble of spins as a qubit. 

5. Nitrogen vacancy centers in diamond 

Figure[T]j): This case is similar to the electron spin just 
discussed, except that the system is a spin triplet with 
5=1, and the ms = state is shifted relative to the 
tos = ±1 states due to crystal-field and magnetic-dipolc 
interactions by 2.88 GHz (59|. The axis of this internal 
field is fixed by the crystal, so that the alignment of the 
states cannot follow the external drive field direction for 
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weak fields. However, if we align the resonator mode 
field along the NV center's crystal axis, and use a drive 
field rotation plane that also contains this axis [fig, (He)], 
we can still realize the desired effect (we take as above 
B c i = 1000G). When the drive field Bd is along the crystal 
axis, the resonator field produces a linear Zccman shift of 
the ms = ±1 states; when the drive field is perpendicular 
to the axis, it mixes all three sublevels, producing a large 
enough avoided crossing that the ms = ±1 states can be 



nearly adiabatically transformed into each other by the 
drive, and we can use them as our two qubit states (Note 
that in experiments, the S = 0, m s = state is typi- 
cally used as \g)). As above, this can be fully accounted 
for using an effective field associated with the rotation. 
As long as the drive frequency is not too large, and the 
drive is as above turned on and off slowly, nonadiabatic 
transitions can be neglected at the error levels of interest 
here. 
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